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Abstract 

We investigate a modification in the action of inflaton due to noncommutativ- 
ity leads to a nonstandard initial vacuum and oscillatory corrections in the initial 
power spectrum. We show that the presence of these oscillations causes a drop 
in the WMAP \ 2 about A% 2 ~ 8.5. As a bonus, from the parameter estimation 
done in this work, we show that the noncommutative parameters can be precisely 
bound to 10 16 GeV or 10 4 GeV depending on the inflation scale. 
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1 Introduction 



It has recently been emphasized that the effects of trans-Planckian physics might be 
observable on cosmological scales in the spectrum of Cosmic Microwave Background 
(CMB) radiation [1]-[T8]. An especially intriguing theory that could makes this possible 
is inflation [19|-[23j. The mechanism of inflation answers several questions that can not 
be solved in the standard big bang cosmology. It is also the first predictive theory for 
the origin of structures in the Universe and these predictions have been verified to great 
accuracy using CMB anisotropy experiments. 

The phase of inflation must have lasted about 60 e-foldings so that the observed 
structures in the Universe had enough time to be seeded from vacuum fluctuations. 
Consequently it can be possible that inflation is begun from scales with characterized 
wavelengths much shorter than the Planck length lp. Therefore one can trace the trans- 
Planckian effects on the inflationary predictions. For instance, due to trans-Planckian 
effects on the vacuum fluctuations, the scale invariant primordial power spectrum Pq 
could be modified to a power spectrum which contains superimposed oscillatory terms 
suppressed by the scale a n = (-f-) n+1 [I21 [2U [251 [26] . Here, H is the Hubble scale during 
inflation and A is the scale associated with the trans-Planckian physics. 

Although the power spectrum and CMB fluctuation can be influenced by the details 
of high energy scale physics, there is no agreement on the size of these effects and n 
can be n = 0, 1, 2 (see [121 I2U [23 [26]). For example using the low energy effective field 
theory and general decoupling arguments, Kaloper et al. [26] have shown a suppression 
of the order of o\ = (-^) 2 and have discussed that it is too small to be observed during 
future CMB experiments. This result strongly violates the first conclusions that the 
deviation would be testable [H [2]. The argument of [26] has been criticized by some 
authors [281 121]. Brandenberger and Martin in |28j have shown that instead of [26J, 
the trans-Planckian physics can leave imprints on the CMB anisotropy. This approach 
is the same as the Danielsson approach [T2]. Both insist on this fact that the vacuum 
state must be modified due to trans-Planckian physics. Although the preliminary data 
analysis [301 ED E2] showed that these oscillatory modifications decrease the WMAP 
X 2 but recently Groeneboom and Elgaroy, by investigating the Danielsson conclusion 
[12] , have discussed that there is no significant evidence for the simulated data to prefer 
the trans-Planckian models [33]. We have also analyzed the Danielsson formula and 
observed no improvement in the \ 2 - 

In this work we change the kinetic part of an inflaton action by assuming a harmonic 
oscillatory term which comes from noncommutativity to resolve the UV/IR mixing 



2 



problem. This term plays the role of a barrier potential and makes the vacuum state 
to be a combination of negative and positive modes. We show that such modification 
of a vacuum state leads to a power spectrum containing oscillatory corrections similar 
to the Danielsson result but with an important difference. Using the cosmological 
Monte Carlo (CosmoMC) code, we will show that due to the presence of this difference, 
our model gives a drop in the y 2 in spite of the Groeneboom and Elgaroy's negative 
conclusion. Having parameters estimated by CosmoMC, the noncommutativity scale 
Ajvc is determined as a function of inflation scale H. Knowing a variety of bounds on 
H, we find the scale of noncommutativity to be of the order of 10 16 GeV or of the order 
of 10 4 GeV depending on what scales for inflation are used. 

The paper is organized as follows. In sections 2 and 3 we briefly review the standard 
calculation of the inflationary power spectrum when the trans-Planckian effects are 
taken into account. In section 4 we derive our modified formula for the power spectrum 
when an extra term is added to the free part of inflaton action due to noncommutativity. 
Comparison of our theoretical predictions with the WMAP data is given in section 5. 
In this section we bound the noncommutative scale using parameters estimated by 
CosmoMC. 

2 Initial Vacuum State 

The authors of [28] have suggested that the evolution of the inflaton modes can be 
separated into three phases I, II and III shown in Figure 1. During phase I, the physical 
wavelength A of modes is smaller than the scale A m lp and the effects of short distance 
physics are expected to be important. During phase II, the physical wavelength is larger 
than the Planck scale but smaller than the Hubble radius H~ l . Since during phase 
II the background is time-dependent, there is a nontrivial mixing between creation 
and annihilation operators at different times[34J. This kind of mixing takes place via 
Bogoliubov transformation which through it a mode v(rj) must be a linear combination 
of positive and negative frequency initial modes [34] 

where a& and f3k are the Bogoliubov coefficients which satisfy the normalization condi- 
tion 

Ki 2 -i^r = i (2) 



3 




Figure 1: Periods 

Finally in region III, modes cross the Hubble radius if -1 and squeeze. The coefficients 
otk and Pk are determined by matching (0Q) to the vacuum state. But in any expanding 
universe such as inflationary period, the notion of a vacuum is ambiguous [31]. The 
reason is that H is time-dependent during inflation and so inflating space-time is not 
exactly a de Sitter space-time in which defining vacuum state is possible. This is a 
well-known problem in curved space-times where the concept of the vacuum state is 
quite ambiguous due to the absence of Killing vector fields [31] . 

A space-time such as de Sitter space-time, admits a timelike Killing vector field, then 
this provides a natural way to distinguish positive and negative frequency modes and 
then similar to the standard procedure in Minkowski space, associate these distinguished 
modes with annihilation and creation operators. By definition, the vacuum state of a 
de Sitter space-time, (known as the Bunch-Davies vacuum [35]), will be the state that 
is annihilated by all the annihilation operators. This vacuum state is invariant under 
the symmetry group of the space-time. For inflationary theories although H is time 
dependent and consequently the Killing vector cannot be defined, this ambiguity can 
be ignored by choosing the adiabatic vacuum when the wavelength of a mode is much 
shorter than the curvature scale of space-time. This vacuum is known as the adiabatic 
vacuum [34J. 

In addition to this, there is also an extra ambiguity as the result of trans-Planckian 
physics below the certain cut-off A through which the vacuum is ambiguous since the 
physics of region I (Planck scales) is unknown. Danielsson approach [12] tries to solve 
this kind of ambiguity. In this approach there is no explicit assumption about the 
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trans-Planckian physics of region I and the emphasize is made only on the point that 
the wavefunctions of the fluctuation modes are not in their vacuum states when they 
enter into phase II [12] . Then for avoiding the ambiguity of region I which is the region of 
unknown physics, one defines the vacuum at the time r/o when the physical momentum 
p of a mode equals the scale of new physics i.e. A [12]. The physical momentum p and 
the comoving momentum k are related through 

k = ap = (3) 
By imposing the initial conditions when p = A, we find the conformal time 770 as 

Without knowledge of the physics beyond the scale A, in order to choose the vacuum 
it is not necessary to take the limit t]q — > 00 but instead one had to stop at the value 
of conformal time given by (j3J). So the question about the trans-Planckian physics is 
that of choosing the vacuum state which clearly would be different from the adiabatic 
or Bunch-Davies vacuum. An alternative vacuum leads to the corrections of order 
a n = (-fO n+1 in the power spectrum. 

3 Region II 

In this section we recall what happens in region II during the inflation. During inflation 
the action of the massless scalar inflaton field tp minimally coupled to gravity is written 
as OS- [23] 

T 

where g^ v is metric and V((p) is the inflaton potential. Usually one assumes a spatially 
flat, homogeneous and isotropic background with the metric 



S = I d 4 x^g{ -g^ipd^ - V(<p) ) (5) 



ds 2 = a 2 {rf) (drf — dx • dx) (6) 

where 77 is the conformal time and a (77) is the scale factor. Here we restrict our attention 
to the scalar perturbations in metric which will be gauge dependent. Hence the line 
element in the longitudinal gauge, takes the form 

ds 2 = a 2 (r]) [(1 + 2$(r7, x))dr/ 2 - (1 - 2^(r/, x))dx ■ dx] (7) 

where the space-time dependent functions $ and $ are the two physical metric degrees 
of freedom which describe the scalar metric fluctuations. In the absence of anisotropic 
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stress, the two metric perturbation $ and ^ coincide, $ = ^. During inflation the 
inflaton field can be separated into a classical background and a fluctuating part 



<p(r),x) = <p {T]) + 6(p(ri,x) (8) 

Because of the Einstein equation, the metric fluctuation $ is determined by the matter 
fluctuation Sip. Thus every thing can be reduced to the study of a single gauge-invariant 
variable field [251 EE] 



a(5 V +^)=zU, (9) 



n 

where a' = da/drj, TC = ~, z = and 1Z denote the curvature perturbation in 
comoving gauge J25J [28]. It is convenient to work with the variable fi s defined by 
fj, s = —\/2kv where k = 8n/mp. After quantizing the theory, in the Schrodinger 
picture, we can write /i s (?7,x) as 

fM s (r], x) = * / d*k \ UFirj) d k e ik - x + tf£i*fa) at fc e"** 1 , (10) 

(27T) 2 7 L J 

where dk and d)_ k are respectively the creation and annihilation operators. The evolution 
of U^irj) is given by the Klein-Gordon equation of motion in Fourier space 

ur^+f^ 2 - ^jui i { V )=o. (ii) 

In the case of slow-roll inflation limit the quantity — is reduced to [25] 

Z " 2 + 3e (19) 
— , i 12 ) 

z r] z 

in which e is the slow roll parameter. The general solution for Ul (77) is written in the 
form 



Ul z (rj) = ^ M* ( a* <H0 H„(\kn\) + (3 k H* v (\k V \) ) , (13) 

where H u is the Henkel function of the first kind. In the above expression, v is a function 
of e in the slow roll case. The coefficients a k and f3 k are the Bogoliubov coefficients which 
are fixed by the initial condition. Using the asymptotic form of the Hankel functions in 
the limit kr\ ^> 1[37] 

H^tkq > 1) « \\-\kr)\e l ^-^~^\ (14) 

V 7T 

one can find the leading behavior of (jj) 

UlUn) « 4^ e-^-^ + -^L e ik ^~^ +■■■. (15) 
V2k V2k 



6 



Naively, one can choose a boundary initial condition by simply stating that the 
modes in the limit 77 — > — oo are positive frequency modes of the kind of Bunch-Davis 
vacuum which leads to a = 1 and (3 = 0. But as have been widely discussed recently 
[DOE], [7]- [18], by the trans-Planckian considerations, the vacuum state is not a Bunch- 
Davis one. To solve this ambiguity, Danielsson [12J assumed the following boundary 
condition 

? K/a) (770) = -ik (U?/a) (vo) (16) 



d 

dr) 



which gives the coefficients a and (3 and consequently the power spectrum P(k) as 



P(k) = P (k) 



1+ < 



sin 



2 ( k 



(17) 



where £ = ^ and 



N . K is a phase factor. A similar result reported in [38J and also see 
and [10] for another discussion and conclusion. In the following section we will 
drive such a formula for trans-Planckian power spectrum by considering an assumption 
about the physics in the region / instead of boundary condition (fl6|) . Our result for the 
power spectrum will be different from (TIT)) and the difference is that in our case two 
coefficients £1 and £ 2 appear instead of £. 



4 Region I 

Since in region I the complete theory of quantum gravity is still unknown, one way to 
tackle with it is to work with effective scenarios. For instance, it is possible to include the 
higher dimensional terms to the Lagrangian [261 [29], [36] or invoke the modified dispersion 
relation models [28 J • Another possibility that has been greatly studied in recent years is 
the noncommutativity of space-time [3] . Noncommutative space-time emerges in string 
theory configuration and is an effective picture of the foamy space-time above the energy 
scale Atvc* (for review see [JHH2]). Field theories on a noncommutative space-time 



[3f,x v ]=i6p, (18) 
are defined by replacing the ordinary product of fields by the star product [131 S3] 

0i * M*) = J ^^0i(^)0 2 (g)e^^^e^ fe+ ^ , (19) 

where the 6^ is a real antisymmetric tensor which introduce a scale K^c = 77^ ■ In Eq. 
(j!9p the tilde denotes the Fourier transform. Because of this star-product, a puzzling 
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mixing of UV and IR scales known as UV/IR mixing appeared during performing loop 
calculations in noncommutative field theories [IS]. The UV/IR mixing makes theories 
nonrenormalizable. Several attempts have been done in order to solve this problem 
|4"6l |4"TI |4"8] . For the case of scalar theories, Grosse and Wulkenhaar in [16], have shown 
that one can cancel out the UV/IR mixing problem by introducing a modification in 
the free part of the action. They assume a harmonic oscillatory term to the action such 
that the kinetic term becomes 



Skin[<t>\ 



d 4 x -d, 



O 2 

+ y5Y 



(20) 



20-V. 



Furthermore in (1201). for the case of Q 2 = 1 the 



where Q G ( 0, 1] and x A 
action is expected to enjoy a duality transformation between position and momentum 



</>(p) «-> vrVl det0| <j>{x) 



x 



(21) 



where 4>(p a ) = J cfix e^ -1 )" 1 ^^ <fi(x a ) . 

Inspired by this modification, we suggest that the kinetic term of an inflaton action 
must be modified due to quantum properties of space-time on the scales where trans- 
Planckian effects become important. Since the noncommutativity is one of effective 
pictures of the trans-Planckian physics, hence for the inflatons, a modified action such as 
(|20l) is considered. Therefor we consider the harmonic oscillator term ^Y(p{x)x 2 (p(x) as a 
correction to the kinetic part of the inflaton action (jSJ), in region I. The noncommutative 
parameter 9^ u is also restricted to the following special case where as the result of it 
the equation of motion of the inflaton is simplified 



e 



( 6 \ 


0-000 

V -6 / 

Consequently under this assumption, x 2 is given by 



(22) 



x 2 



(2^>1 2 = 4^ [( 



x 3 ) 2 - (x 1 ) 2 



(X 



2\2 



- m 



and the kinetic part of the inflaton action is written as 

S km = j d'x^ Q sTd^ d uV + 2^ [(x 3 ) 2 - (x 1 ) 2 - (x 2 ) 2 - (r?) 2 ] V 2 



(23) 



(24) 
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Now the variable n s (rj,x) is decomposed into the Fourier modes in the following 
general form 

fx s ( V , x) = * f [dk] ( a k U k (rj) V„(x) + a\ U* k V?(x) ) , (25) 

where [dk] is the measure of the integral and the V k (x) is the eigenfunction of the three 
dimensional Laplace operator A3 

(A 3 + (x 3 ) 2 - (x 2 ) 2 - (x 1 ) 2 ) V k (x) = di (V^gpdj) 14(x) 

+ ((x 3 ) 2 - (x 2 ) 2 - (x 1 ) 2 ) V k {x) 
= fcV fc (x) , (26) 

in which the metric g % i has been defined in ([6]) and k 2 is the eigenvalue. Also in the 
expansion fl25l) . Uk(r)) obeys the following equation 



a" 



U h {rt) = 0. (27) 



In the case of the absence of the second term in equation ( |26l) . the 14 (x) are simply 
plane- wave e lk x and the integration measure in (125]) for this case is f[dk] = J d 3 k. The 
differential equation (126]) can be solved by separating Vjt(x) into 

Vfc(x) = V r fe3 (x 3 )Vfc 1 , fc2 (x 1 ,x 2 ) (28) 

The differential equation for the Vfc 3 (x3) reduces to the harmonic oscillator differential 
equation, which its solution is given in terms of Hermite polynomials. Thus, the wave 
number ks takes discrete values which we label by positive integers k > 1. Also for the 
14i,fc 2 (xi5 X2) we find a differential equation, which its solution is known as parabolic 
cylinder functions [37J. The parabolic cylinder function, D p (z) with continues p, is a 
class of special functions defined as the solution to the following differential equation 

^^+(z 2 + X( P ))D p (z)=0, (29) 
together with the integral representation as 

DJ Z ) = — / e- xz ~x- p - l dx (Rep<0). (30) 

r(-p) Jo 

Since k 3 is discrete and k\ and k 2 are continuous, the integration measure in (|25j) is 
given by 

POO POO 00 

[dk] = / d 3 hd 3 k 2 V . (31) 
Jo Jo fe3=1 
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Because the modes are in region I, one expects — — > for the time-dependent part of 
modes in ( 1271) . Hence equation ( 1271) is simplified to 

K{t) + (k 2 + r 2 )U K (r) = , (32) 

in which Uk{fj), fj and k have been redefined as follows, 

U k (t) = ( ^ J U k ( V ), r = k 2 = V29k 2 . (33) 



The time evolution of modes in ( 1321) is equivalent to the following situation. Consider 
a scalar inflaton on a non-static background 

ds 2 = C(rj) (dr] 2 - rfx ■ dx) . (34) 

It can be shown that with a conformal scale factor representing a bouncing universe 

C{r[) = a 2 + b 2 r] 2 - oo < 7] < oo, (35) 



one obtains the same equation for the evolution of time-dependent modes [Mj. The 
exact solution to the equation (1321) is then given in terms of parabolic cylinder functions 



Ui{r) = U K (r) = D 1+iK * ( (1 + i)r ) . (36) 

2 

Now we study the asymptotic behavior of U^(r) = U R {r) modes living in region I 
in order to match them with U^(t) modes living in region II at the boundary r] . This 
helps us to obtain coefficients a>k and [3k- The parabolic cylinder function D p (z) has the 
following linear relation [37J 

D p (z) = e-*** D p (-z) + ^— - e--^~ D^(iz) , (37) 

where z = (1 + i)r and p = — x+% jf ■ This relation is a useful relation in order to rewrite 
the modes U^(t) as a combination of positive and negative modes. The asymptotic 
behavior of D p (z) for the large values of p is 



2 4+2 / f\ p\ f 1 p\ I / 1 p\ z 2 



(8("1 - I)) 



DM* : . -TT exp|| ; + | h -2. - i + ^ - |. (3* 



10 



Substituting fl38|) into the relation fl37|) with p = — l K , we find 



-P i+j^((l+i)r) 



2"^ e 



+ ™ Uri_t« In ■- 



2e ^ k 2 



jlt 2 / 1 in \ -7Ti 

2tt 2~ eV 2 4 J 2 



2\ „„- 2 -2 



r (§ + <£) 



2e 2 k 2 



(2k 2 ) i 



27T fiT7T 



ex P 



ZV27T 



r(| + if) (2« 2 )i 



exp 



8 4 



(39) 



Thus Ulirf) can be written as 



— a k e- ik « + -^=(3 k e ik \ 



(40) 



/2k V2k 

where the coefficients a k and (3 k are found by matching the solution fl39|) with U^ij]) 
in ffl~5l) at ^ = ^ as 

ITT k^7v9 

a k = ie 8 2v/2 e * 

and 

P k = 



Using the identify [37] 



1 1 2 + * 2 >/2 



I I 2 + % V 



(41) 
(42) 



7T 



cosh Try 



one can verify the normalization relation 



\a\ 2 -\3\ 2 = l 



In general the gamma function can be represented in the following form 



i ( l - + ik 2 e 



\ ( - + tk 2 e 



71 



cosh 1/2 7rk 2 6 



(43) 



(44) 



(45) 



where (f(k 2 8) is a phase factor. We will employ this relation for the next calculations 
of P(k). 
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Using equations (J41j) and (J42|) . one can calculate the corresponding energy density 
and pressure by the mean value of energy- momentum tensor (T^) in the trans-Planckian 
region [27] 

2 

dkk 2 I a 2 



(P) 
<P> 



1 



4vr 2 a 4 
1 

4vr 2 a 4 



d&A; I a 




(46) 
(47) 



Inserting the mode function Ul(r]) into the vacuum expressions of the energy density 
and pressure and using the normalization condition (j4"4"l) . one finds 



(P) 



47r 2 a 4 



^ 3 |/3 fc | 



<P> 



47r 2 a 



(48) 
(49) 



in the trans-Planckian region. Now at this stage we turn to the superhorizon scales 
{krj << 1) and start to compute the power spectrum. First we substitute ( 14T1) and ( 1421 ) 
in (Tl3|) and take the superhorizon limit. Since in this limit the Hankel functions take 
the following form 



H v (kr) < I) 



2 ,vr „ 3 r(z/) ., , 



(50) 



The amplitude of f/^ 7 modes is given by 



117, 



ii i 



2"-f i» 
2lr(|) 



\oih + Pk\(krjY 



(51) 



which gives the scalar power spectrum as 



P(k) 



k 3 



2tt 2 



ul 1 



I +2/? ( ( — 



I/3P 



a: 



= Po(*0 K| 2 

where Pq is the ordinary power spectrum for the scalar modes given by 
Usually Pq is parameterized by 



(52) 



(hV [ k - 3 - 2u 



\2n J \aH 



(53) 



1l a -l 



(54) 
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in which n s is the conventional definition of spectral index and A s is the spectral ampli- 
tude for the scalar perturbation, k* is a scale which is fixed to be 0.05 Mpc~ l [50]. For 
the case of slow roll inflation, the Hubble parameter H is A;-dependent with the form 

m 

H ( k x ~ 



H \h 

where e is slow-roll parameter. Then the quantity krj varies as 



(55) 



A NC Anc ( kY 

k " = ^r~m{Tj ■ (56) 

Similarly, the quantity k 2 9 varies as 

. -i 2 

(57) 



,2^2 

ee ' 



,20 a P 



A 



NC 



A NC ( k x 



Hn \k 



We now substitute gTfl and (|32jl in (|52|) and use (|5Tj) and (Pj) to find P{k) as 
P(k) = A s 

e~ + V2e ^ cosh 1/2 j^^=r ) cos(vr/4 + 2A;r/) + l + e ~)(58) 



where 



e -c*-f = e -^-B5-; ^ . (59) 

Since the /c 2 # is large in value, the exponential terms decay very fast. Thus it is possible 
to approximate the exponentials with a polynomial function as follows 



-ck 2 6 -c 



'^nc\ 2 ( k \~ 2s 



Keeping only terms up to order and redefining the coefficient cij^ as a new 
parameter £i we will find 

p ^ - a - (ST x + 6 ( £ r cos ( V4 + 1 ( ^ r ) + - ■) • (6i) 

where £2 = ctq = a^t- This result that the coefficients £1 and £2 in the power spectrum 
formula are not equal, is the difference between our result and Danielsson formula and 
will play a crucial role in CMB data analysis. 
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We saw that by considering noncommutativity corrections to the free part of the 
action of inflatons in region /, the power spectrum receives oscillation corrections. This 
result has a simple physical interpretation. With a suitable changing of variables, 
equation f|32|) for the primordial modes is similar to the Schrodinger equation for a wave 
function ip(x) of a quantum mechanical particle in a one dimensional barrier potential 
V(x) = x 2 

^ + (E- V(x))il> = 0. (62) 

Because of this kind of barrier potential we expect modes in region I be a superposition 
of the incoming modes and the reflected modes which have been scattered off the barrier. 



5 Effects on CMB Temperature Fluctuations 

In this section a CMB data analysis is provided in order to show that the agreement 
of our modified power spectrum with the CMB data is better than the standard one. 
The fluctuations in the temperature of CMB radiation can be expanded in spherical 
harmonics pj2, [20] 



AT(d,<Pj ^ 



oo m=l 

^ ai m Yi m (6,4>) (63) 



rp 2—t 

1=2 m=-l 

for convenience, we have excluded the monopole and dipole terms. The initial power 
spectrum P(k) given in (|6T|) is related to the CMB anisotropy through the angular 
power spectrum G\ which is defined by two-point correlation function of the temperature 
fluctuation 

AT . A . AT . A . \ ^ 2/ + 1 



m — {n 2/ 



T T ^ 4vr 

1=2 



-^-C, PUcobO), (64) 



where Pi(cos8) is the Legendre polynomials and cos 9 = h\ -n^. The fi\ and fi2 are unit 
vectors pointing to arbitrary direction on the sky. The angular power spectrum C\ can 
be related to P(k) through 

/°° dk 
T 2 (k)P{k) — (65) 

in which T(k) is the transfer function. Since Thomson scattering polarizes light [52J, 
there are also angular power spectrum coming from the polarization. The polarization 
can be divided into a curl (B) and curl-free (E) component which yields four independent 
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angular power spectrums as Cf T \Cf E ', Cf B and the T — E cross correlation Cf E [52J. 
The WMAP experiment has reported data only on Cj T and Cf E [52J. The previous 
data analysis has reported different conclusions about the signatures of trans-Planckian 
modification of power spectrum [30j-[33j and [53]- [59]. The analysis of [3QI [311 [32] 
shows that the presence of oscillations in power spectrum causes an important drop 
in the WMAP x 2 °f about A% 2 « 10. \ 2 is defined as \ 2 = ~ 2\og C(x\p) in which 
C(x\p) denotes the likelihood which is a conditional probability function that allows 
us to estimate unknown cosmological parameters p based on CMB data x and satisfies 
the normalization condition J C(x\p)dx = 1. In spite of the conclusion of [301 13TI I32J. 
the recent analysis of [33] claims no evidence for trans-Planckian oscillations especially 
in the Danielsson model. Now we are going to verify our prediction for the trans- 
Planckian power spectrum by analyzing the WMAP data. In the previous section we 
found trans-Planckian oscillatory corrections to the power spectrum (|6T!) in which the 
coefficients £i and £2 are not equal in comparison to the Danielsson formula [12J. In 
order to compare our result ( ISTT) with the recent WMAP data, we present a Bayesian 
model selection analysis data using the CosmoMC (Cosmological Monte Carlo) code 
developed in [50] which makes use of the CAMB program [60]. This program employs 
a Markov-Chain Monte Carlo (MCMC) sampling procedure to explore the posterior 
distribution. A Bayesian analysis provides a coherent approach to estimating the values 
of the parameters, p, and their errors and a method for determining which model, M, 
best describes the data x. Bayes theorem states that 

PW,M) ^ p ^ M \ (66) 

where P(p\x,M) is the posterior, P(x\p,M) is the likelihood, P(p\M) is the prior and 
P(x\M) is the Bayesian evidence. Conventionally, the result of a Bayesian parameter 
estimation is the posterior probability distribution given by the product of the likelihood 
and prior. The parameter space we consider is h (the dimensionless Hubble constant), 
flj, (the amount of baryons), Q c (the amount of cold dark matter), (the amount 
of dark energy), r (the redshift of reionization), n s , A s , £1, £ 2 and e. We altered the 
CosmoMC to include the new parameters £1, £2 and e. To find the best-fit values for the 
parameters, we use the recent five-year WMAP (WMAP5) dataset [61] . We find that our 
model gives a better fit with A% 2 ~ 8.5 as compared to the standard inflationary model 
without trans-Planckian corrections. This improvement in the value of WMAP \ 2 has 
been obtained because of the difference in the £1 and £2 coefficients. We also checked 
that there is no improvement in \ 2 in the Danielsson model with equal coefficients 
£1 = £2- In Fig. (T51) we have plotted the ID marginalized posterior probability and 
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A s x 10 10 
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X 2 /d.o.f. 







0.71 


0.023 


0.11 


0.73 


0.083 


21.32 




0.96 


2658.48/1453 


0.309 


8.14 


0.72 


0.022 


0.11 


0.75 


0.089 


21.03 


0.0798 


0.95 


2649.96/1456 



Table 1: Best fit parameters from the WMAP data for the standard model of inflation 
compared to the best fit parameters obtained for our trans-Planckian model. 



the normalized mean likelihood for all of the primordial parameters of our model. In 
Fig. (131) the contours of the 2D marginalized posterior probability and the normalized 
mean likelihoods are displayed for each parameters pairs. The shading shows the mean 
likelihood of the samples and helps to demonstrate when the marginalized probability 
is enhanced by a longer parameter space rather than by a better fit to the data. Using 
CAMB |60j we have plotted the TT and TE angular power spectrum in Fig. (j4j) and 
the matter power spectrum in Fig. ([5]), corresponding to our trans-Planckian model. 
The presence of small oscillations in the trans-Planckian angular TT power spectrum 
for / < 200, is the reason that we obtain a smaller value for the likelihood \ 2 ■ 

Table CD shows the best fit parameters found by CosmoMC for the case £i = (stan- 
dard inflationary model) and ^ ^ (trans-Planckian model). For the case of £i = 
the standard values for the cosmological parameters [62] are recovered. We can esti- 
mate the scale of noncommutativity, ^= = Ajvcs from the trans-Planckian parameter 
£ 2 = H/Anc ~ 10~ 4 . If the scale of H during inflation is approximated to be of the 
order of H > 10 12 GeV [63], then the noncommutative scale is estimated as A^c > 10 16 
GeV. On the other hand as Refs. [Ml EHJ EHl [67] have demonstrated, the Minimal 
Supersymmetric Standard Model (MSSM) has all the ingredients to give rise to suc- 
cessful inflation. The MSSM inflation occurs at low scales i.e. H > 1 GeV. Then this 
put a bound on the scale of noncommutativity as A^c > 10 4 GeV. Similar bounds on 
noncommutative scale, has reported in [6H] using the investigation of Lorentz symmetry 
violation due to noncommutativity during experiments and also in [69J by studying of 
hydrogen atom spectrum and the Lamb-shift effect on a noncommutative space which 
is A NC > 10 TeV. 

6 Conclusion 

We showed that a modification in the action of one inflaton due to noncommutativity, 
leads to a nonzero Bogoliubov coefficient (5k which affects the power spectrum to have 
oscillatory corrections. Using CosmoMC code we showed that these corrections cause a 
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Figure 2: One dimensional marginalized posterior probability distribution (solid lines) 
and mean likelihood (dotted lines) for the cosmological parameters. 
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Figure 3: Posterior constraint for our model. The contours show the 68 % and 95 % 
confidence limits from the marginalized distribution. 
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Figure 4: The best fit obtained for our trans-Planckian angular TT (up) and TE (down) 
power spectrum compared to the best fit standard inflationary power spectrum and five 
years WMAP data. 
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Figure 5: Predicted matter power spectrum for our trans-Planckian model compared 
to the 2dF data [70]. 
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drop in the WMAP x 2 of about A% 2 ~ 8.5. We also derived the noncommutative scale 
using parameters estimated by CosmoMC. Tacking the scale of inflation to be H > 10 12 
GeV, we can estimate the scale of noncommutativity which yields A^c > 10 16 GeV and 
with the scale of inflation H > 1 GeV, we get A^c > 10 4 GeV. 
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